Elementary excitations in one-dimensional spin-orbital models: 
neutral and charged solitons and their bound states 
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We study, both numerically and variationally, the interplay between different types of elementary 
excitations in the model of a spin chain with anisotropic spin-orbit coupling, in the vicinity of the 
"dimer line" with an exactly known dimerized ground state. Our variational treatment is found 
to be in a qualitative agreement with the exact diagonalization results. Soliton pairs are shown to 
be the lowest excitations only in a very narrow region of the phase diagram near the dimer line, 
and the phase transitions are always governed by magnon-type excitations which can be viewed as 
soliton-antisoliton bound states. It is shown that when the anisotropy exceeds certain critical value, 
a new phase boundary appears. In the doped model on the dimer line, the exact elementary charge 
excitation is shown to be a hole bound to a soliton. Bound states of those "charged solitons" are 
studied; exact solutions for JV-hole bound states are presented. 
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I. INTRODUCTION 

One-dimensional systems with coupled spin and or- 
bital degrees of freedom have been attracting a consid- 
erable attention during the last two years. This interest 
is partly motivated by the progress in the experimen- 
tal study of the quasi-one-dimensional spin-gap mate- 
rials Na2Ti2Sb20 and NaV205, which are believed to 
be described by a two-band [Ocbitally degenerate Hub- 
bard model at quarter filling Im Orbital degrees of free- 
dom may be described with the help of pscudospin-1/2 
variables; the corresponding effective Hamiltonian for the 
two-band Hubbard model was derived long ago by Kugel 
and Khomskiicl and has the following form: 

71 

+ £ {Js{S n S n+1 ) + J T {r n T n+l + £« +1 )} , (1) 

n 

where S n and t„ are respectively the spin-^ and 
pseudospin-i operators at the site n. The Hamiltonian 
(|l|) can be viewed as two spin-^ chains coupled by a bi- 
quadratic interaction only (a generalized spin ladder). 

Generally, the model has an SU(2) x U(l) symmetry, so 
that the total spin (S, S z ) of the S-chain, as well as the 
z-projection r z of the total spin of the r-chain are good 
quantum numbers. Under certain simplifying assump- 
tions (neglecting Hund's rule coupling, nearest neighbor- 
hopping between the same type of orbitals only, and only 
one Coulomb on-site repulsion constant) , one obtains the 
above Hamiltonian with e = e' = 0, J s ^t-J T = jK which 
possesses the hidden SU(4) symmetry.Q'B At this special 



point, the model is Bethe ansatz solvable^ and gapless. 
Depending on the microscopic details of the interaction, 
this high symmetry can be broken in several ways. For 
example, finite Hund's rule coupling or existence of more 
than one Coulomb repulsion constant makes J s and J T 
different, and local crystal fields can induce considerable 
anisotropy in the orbital sector. 

The phase diagram of the isotropic (e = e' = 0) 
SU(2) x SU(2) version of the model (|l|) in the vicinity of 
the ST/ (4) point was analytically studied recently by sev- 
eral authorsJ3~la Extensive density matrix renorr 
tion group (DMRG) studies of the phase diagramlHH: 
established the form of phase boundaries aruLthe quan- 
tum numbers of the lowest-lying excitations.liil 

Moving away off the ST/ (4) point towards larger J s , J r , 
one runs into the spontaneously dimerized phase with 
a finite gap and twofold degenerate ground state. The 
structure of the elementary excitation spectrum in the 
dimerized phase is known rather incompletely. The weak 
coupling (J S , T 3* K) region of the dimerized phase is 
known to be. a realization of the so-called non-Haldane 
spin liquidE3 where magnons become incoherent exci- 
tations since they are unstable against the decay into 
soliton-antisoliton pairs. Variational study for a spe- 
cial "dimer point" J s = J T = jK, where the exact 
ground state in— a form of a simple valence bond solid 
state is known,ll3 confirms that solitons are the lowest 
excited states also in |Gase of strong coupling. However, 
recent DMRG resultsliil revealed a complicated pattern 
of the change in quantum numbers of the lowest-lying 
excited states in the vicinity of this special point; the 
standard DMRG technique provides only the informa- 
tion about the quantum numbers and energy of the low- 
est excited state, but not the information on the entire 



1 



spectrum. Recent workli-J using the so-called dynamical 
DMRG method contains data for the dynamical struc- 
ture factor, but for only one point in the dimerized phase, 
namely for the dimer point. 

Further, realistic material-relevant models are gener- 
icalWi-anisotropic (e, e 1 ^ 0) and asymmetric (J s 
J r );B'ta the anisotropic model of the same type., arises 
also as an effective Hamiltonian for a spin tubeJia Thus 
we think it is of interest to get insight into the properties 
of the elementary excitations in the dimerized phase of 
the more general asymmetric/anisotropic model, which 
is the aim of the present paper. 

In section |n| we introduce briefly the "dimer line" in the 
space of model parameters where the exact twofold de- 
generate spontaneously dimerized ground state is known. 
Based on this picture of the ground state, in section 
III we present the variational study of the interplay be- 
tween different types of elementary excitations, namely, 
between soliton-antisoliton pairs and their bound states 
(magnons). We show that soliton pairs are the lowest 
excitations only in a narrow region of the phase diagram, 
and away from the dimer line the low-energy physics is 
determined by magnon-type localized excitations which 
are essentially soliton-antisoliton bound states. An ad- 
ditional phase boundary is shown to appear when the 
anisotropy exceeds a certain critical value. Variational 
dispersion curves are found to agree qualitatively with 
the exact diagonalization (Lanczos) results in the vicin- 
ity of the dimer line. 

In section |iv| we describe some exact results for the 
doped spin-orbital model on the "dimer line". We 
show that the natural charge excitation in this case is 
a "charged soliton", i.e., the bound state of a hole and 
a soliton in the dimer order. It is found that under cer- 
tain condition "charged solitons" themselves exhibit a 
tendency to form bound states; exact expressions for the 
energy of those bound states are presented. 



II. EXACT GROUND STATES 



where we have set K to be the unit energy scale. Then 
it is easy to see that the choice J T — 3/4, J s = (3 + e)/4 
makes the checkerboard-type singlet product wave func- 
tion an eigenstate, and one can prove that it stays the 
ground state as long as —2 < e < oo. Thus, in the 
3D phase space ( J s , J Tl e) one has a line with an exactly 
known ground state; in what follows we will call it the 
"dimer line" for the sake of brevity. This line ends at the 
multicritical point e — — 2, where the ground state be- 
comes degenerate with infinitely many other states: e.g., 
changing any of the orbital (t) singlets into a triplet 
with t z — ±1 does not change the energy. The point 
e = +oo is another multicritical point where the triplets 
with t z — can be created with no energy i-cpst. For 
e < — 2 the ground state is also exactly knownEZI and is a 
product of the ferromagnetic state in r-sector and the an- 
tiferromagnetic Bethe-ansatz state in the S-sector. The 
energy of ferromagnons softens at k = tt when e tends to 
-2. 

We have also checked that (^) at J T = 3/4, J s = 
(3 + e) /4 is the only nontrivial model with exactly known 
ground state (within the physically relevant subspace of 
Hamiltonians defined by (|lj)) which can be found within 
the matrix-product based optimum ground states ap- 
proach (see e.g. Ref. [l| and references therein), at least 
if one uses a rather general matrix-product ansatz 

|*) = Tr( 9l (C)g 2 (C) ■ ■ ■ g N ^{C)g N (C)) 

where C = {x,y,x',y',p,q}, and the elementary matrix 
g n {C) has the following form: 



9n(C) 



X\ U)r 
q\ TT>r, 



y\ IT)r 



-P\ U)n 

y'\ 11% 



Here n denotes the rung of the "ladder," and the first 
and second arrows denote the states of S and r spins, 
respectively. 



In addition to the SU(4) model, there are several 
other particular points in the phasg diagram of ([!]), for 
which exact results are available .EaEZliiS Particularly, at 
the point J s = J T = 3K/ 4 the exact ground stated is a 
product of checkerboard-ordered spin and orbital singlets 
(see Fig. [|) , so that the spin and orbital degrees of free- 
dom are completely decoupled in the ground state. This 
ground state is spontaneously dimerized and thus twofold 
degenerate in case of periodic boundary conditions (for 
open boundarieSpthe degeneracy would be 8-fold). 

One can showtH that the same ground state persists 
in a finite range of anisotropy in the case of the special 
choice £ = £■'. In this case, up to a constant term, the 
model ([l]) can be recast in the form 

H = Y,(SnS n +l + Jr){T n T n+1 + + J s ) , (2) 



III. VARIATIONAL STUDY OF THE 
ELEMENTARY EXCITATIONS 



In the vicinity of the dimer line, the elementary excita- 
tions can be studied variationally. There are essentially 
two types of excitations: (i) solitons connecting the two 
degenerate spontaneously dimerized ground states (ex- 
actly speaking, soliton pairs), and (ii) magnons which 
are localized excitations and can be viewed as soliton- 
antisoliton bound states. The variational results give an 
upper bound for the dispersion at the dimer line where 
the ground state is exact, and one may expect that they 
provide a reasonably good approximation in the vicinity 
of this line. 
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A. Solitons 



For a single soliton excitations one can use the varia- 
tional ansatz of the type introduced in Rcf. [l3| 

= E{ eipX "l L >« T +ae*^\R)^} . (3) 



Here the states \L)° T , \R)n T are shown schematically in 
Fig. 0, and x^ = 2n + |, x„ = 2n + | are the effective 
"center of mass" coordinates of \L) and \R) states, re- 
spectively, and a is a variational parameter. This soliton 
state has the quantum numbers S z — er, r 2 = t and may 
be combined in pair with the corresponding antisoliton 
state to get states with (S,t z ) = (1,±1), (0,0), (1,0) or 
(0,±1). Since there are no diagonal exchange interac- 
tions in our "ladder" , the variational energy of a single 
soliton does not depend on (er, r) so that we omit the 
corresponding index in what follows. 

It is easy to see that there exist a couple of useful sym- 
metry properties. Upon the operation of interchanging 
upper and lower leg (S «-> r) the R and L soliton states 
transform as \L) n h-> \R) n -i/2, \R)n >-> \L) n +i/2, so that 



\a;p) i-> a\(a) 1 ;p) 



(4) 



That means that if the variational energy of the state (||) 
calculated with the Hamiltonian (||) has a minimum at 
a certain a — ao, then the energy of the same state cal- 
culated with the Hamiltonian with interchanged S and r 
operators should have a minimum at a' — 1/a. An im- 
portant consequence is that in the symmetric case e = 0, 
J s = J T the minimum should always occur at a = ±1. (It 
should be remarked, that this property was missed in our 
earlier treatment ;£3 the variational parameter £ = ae zp 
used in Ref. [l3|was erroneously assumed to be real, which 
forced the minimum to be at £ = ±1). 

Upon a reflection x t— > — x the soliton states transform 
as \L) n i— > |i?)jv_ n , \R) n i— > \L)pf- n , so that one has 



\a;p) 



a^ 2N+1 )\{a*)- x ;p)* 



(5) 



where the asterisk denotes complex conjugation. The en- 
ergies of two states with mutually conjugate wave func- 
tions should be equal, therefore replacing a by 1/a* 
should not change the variational energy, which means 
that the minimum is reached on a subspace with \a\ 2 = 1. 
Thus, we can set a — e itp in our variational ansatz (^|). 

Essentially simple but tedious calculation of the vari- 
ational dispersion yields 

F i \ _ Fo (p) + F c (p) cos <p + F s (p) sin <p 
solW 24(25- 16 cos 2 p) (5 + 4 cos y> cos p) ' [ > 

where we have used the notation 



F (p) = [J T (e + 3) + 3J S ](383 - 392cos(2p)) 

+ (e + 3) (48 cos(4p) + 180 cos(2p) - 93) , 
F c (p) = 4cosp(33 - 32cos2p){(2J r - 3)(e + 3) + 6J S }, 
F s (p) = 20sinp(25- 16cos 2 p){J T (e + 3) - 3J S } . (7) 

It is now easy to minimize the expression (j^j in ip and to 
obtain the final variational soliton energy. The equation 
dE s /dip = always has two roots, (po and ir — tp$, and 
one of them yields the minimum for < p < 5, and the 
other does so for ^ < p < it. In this way, though one 
cannot see that from Eq. (^|) directly, the final minimized 
soliton energy £™[ n (p) is symmetric with respect to the 
shift p i — > p -|- 7r, in full agreement with the fact that the 
Brillouin zone in the dimerized phase is halved. 

For certain special cases the variational soliton disper- 
sion can be sufficiently simplified. It is easy to see that 
under the following condition of "generalized symmetry": 



J s = ./(l + e/3), 



J T = J 



(8) 



the function F s (p) identically vanishes, so that the min- 
imum of is reached at ip — or ir (i.e., a — ±1), and 
one obtains 



£ s T = (3 + e)/24 



62J + 12cos2p + 3- (64J+ 12)|cosp| 
5 — 4| cosp| 



(9) 



At the dimer line, i.e., for J 
simplifies further to 



|, the above expression 



rndimcr 
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-(5 — 4| cosp[ 



(10) 



The energy gap at the dimer line occurs always at p = 0, 
7r. One can check that the gap which follows from ( |l0| ) 
coincides with the result given in Ref. |l3|; however, for 
general p the formula ( |To| ) yields the lower energy and is 
in addition much simpler than the corresponding expres- 
sion in Ref. [ll|. 

The energy of the continuum of soliton-antisoliton 
scattering states is given by 

E pair (k, q) = E s (k/2 + q)+ E a (k/2 - q) , 

where k and q are respectively the total and relative mo- 
mentum. All momenta here are defined for the original 
non-dimerized chain, and the Brillouin zone is in fact 
halved, so that k E [— 7r/2, 7r/2]. 



B. Magnons 

We understand "magnons" as localized excitations 
which do not break the dimer order (in contrast to soli- 
tons). In a crude approximation, one could imagine that 
a magnon state corresponds to exciting one of the sin- 
glet bonds on the S or r chains to a triplet and making 
it propagate along the chain. One may speak about S- 
or r-magnons, respectively. In turns out, however, that 
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such a crude single-mode approximation overestimates 
considerably the magnon energy. A simple reason is that 
a state \t ntl ) shown in Fig. || is not an eigenstate of the 
Hamiltonian even at the dimer line, and under the ac- 
tion of the Hamiltonian it gets mixed with other states. 
In the vicinity of the dimer line the largest contribution 
to the first-order matrix elements (i.e., of the {ip\H\t n ^) 
type) comes from the action of the local Hamiltonian h n 
coupling the sites (2n + 1) and (2n + 2) (see Fig. ||). The 
other first-order contributions are proportional to small 
parameters A r = J T — | , \ s = J s — , measuring the de- 
viations from the dimer line, and we neglect them for the 
sake of consistency (see below). Thus, it looks reasonable 
to take a linear combination of the two orthogonal states 
\t niJ ) and 

|A ^n/i) h n \t ri ^j {t n fi\ h n\t n fi) 71 (H) 

(see Fig. ||) as a simplest variational ansatz. The magnon 
state with momentum k can be written as 

|m; k) = Y, e ' lH2n+1) {cosxIW) + sin X |A«t nA1 )} . 

n 

(12) 

where |A^H n/1 ) is assumed to be normalized, and x 1S 
a variational parameter. After minimization in Xi the 
magnon dispersions take the following form: 

E s (k) = e 0iS + (l/8)(3 + e - 4J S ) cos(2fc) , 
E°(k) = e° T + (1/8)(1 + e)(3 - 4J T ) cos(2fc) , (13) 
E±\k) = 4$. + (l/8)(3 - 4J T ) cos(2fc) , 

where the superscripts for r-magnons indicate the value 
of t z . One may observe that the inclusion of the above- 
mentioned contributions to \^^t n ^) which are propor- 
tional to A s , A r would lead to corrections of the order 
0(A 2 T ) in both the bandwidth and the self-energy of 
magnons. To be consistent, we have to omit those con- 
tributions, since we have already implicitly neglected the 
corrections to the ground state energy which are also 
0(Ag T ) in the leading order. (This does not mean, 
of course, that our calculation will have the accuracy 
0(Ag T ), since the contributions from the higher-order 

matrix elements of the type (ip\(H) m t n g) do not contain 
any small parameter. The ansatz (12) should be con- 
sidered as purely variational, sort of "one-step Lanczos" 
trick). 

The self-energies for r-magnons can be compactly writ- 
ten as 

e£ T =3J,/4 + J T 

(l + 7 ? £) + 2{(l + r7£) 2 +48j2} 1/2 
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(14) 



cumbersome due to the fact that the state {Ay-'tn^) in 
case of the S'-magnons does not have a certain value of 
the total orbital momentum T and is a mixture of states 
with T = and 2: 



eo.s — Js + Jt — 



AJ T (3 - e) + 3{1 + e) + VZ 



16(e 2 + 2e + 3) 
Z = 16 J 2 (e 2 + 2e + 2)[e 4 + 4e 3 + 27(e 2 
+ 8J r e 2 (e 4 + 6e 3 + 21e 2 + 32e + 24) 



2e + 2)] 



where r\ — 1 for /i = ±1 and r\ = —1 for fi = 0. The ex- 
pression for the self-energy of an S'-magnon is somewhat 



+ e 6 + 6e 5 + 21e 4 + 44e 3 + 36(2e 2 + 2e + 1) (15) 

One can observe that the ansatz for a magnon state 
which we have used above is in fact equivalent to the 
simplest possible ansatz for a bound state of soliton and 
antisoliton, where only the configurations of the type 

\L) n \L) n+h , \R) n \R) n+i , \L) n \R) n+i , \R) n \L) n+h 

with soliton and antisoliton being nearest neighbors are 
taken into account. Here the bars denote antisoliton 
states; note that for fixed boundary conditions, solitons 
and antisolitons live on different sublattices, e.g. the free 
spin in the r chain (see Fig. |l|) is placed on the odd (even) 
rung for solitons (antisolitons), respectively, which is re- 
flected in the above notation. 

This indicates that magnons should be viewed as 
bound soliton-antisoliton states, rather than as indepen- 
dent excitations. Our simple ansatz can capture only 
the short-range physics of the tightly bound states with 
quantum numbers (S,t z ) = (1,0), (0,0), and (0, ±1); it 
will miss, however, possibly existing higher-lying loosely 
bound states as well as states with (S,t z ) = (1,±1). 

It is interesting to remark that solitons themselves 
can be viewed as bound states of two spinons belong- 
ing to the S and r chains, respectively, which should un- 
bind when one moves towards the exactly solvable SU (4) 
point. Therefore on the way from the SU (4) point to the 
dimer line there should exist an interesting hierarchy of 
bound states, which would be an interesting subject of 
the future studies. 



C. Crossover between solitons and magnons 

The minimum of the lower boundary of the soliton con- 
tinuum is at k = in the immediate vicinity of the dimer 
line. When J s , J T decrease below certain values, the sec- 
ond equivalent minimum splits off and starts to move to- 
wards 7r/2. The magnon gap always lies either at k = tt/2 
or at k = 0. Fig. [| illustrates the behavior of the gap 
along the "generalized symmetric line" J s = J(l + e/3), 
J T = J for a few choices of e. 

One can observe the following important feature: the 
solitons are the lowest excitations only in a limited region 
of the phase space. Even at the dimer line line solitons 
cease to be the lowest excited states if the anisotropy 
is sufficiently strong: the -magnons become the lowest 
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excitations for e < —59/33 and soften at the phase transi- 
tion point e — —2 (respectively, r°-magnons become the 
lowest excitations for e > 59/7 and the corresponding gap 
saturates at the finite value of § at e — > +oo). One can 
see that whenever the variational gap closes while mov- 
ing from the dimer line, the magnon gap touches zero 
first, so that the phase transitions are always governed 
by magnons. 

Three branches of r-magnons are degenerate at e = 0, 
and for s > (e < 0) the r°-magnons (r ± -magnons) 
become lower in energy, respectively; the 5-branch lies 
always in between the and r° branches. In this con- 
nection, we would like to make a remark concerning Fig. 
4 of Ref. which shows the behavior of the S-magnon 
gap (the lowest excitation in the (S z ,t z ) = (1,0) sector) 
across one of the phase boundaries in the model (|J) with 
e = — 1. We think the behavior of S'-gap is not relevant 
as a proof of the first order transition, because another 
excitation, namely that of the (0, 1) sector, should be 
responsible for the transition. 

Away from the dimer line our description should be- 
come progressively worse since the pure dimer state is 
too far from the correct ground state. Thus the present 
approach cannot be used for estimation of the shape of 
the phase boundaries, except when they become close to 
the dimer line (see below). 

Another observation is that though the gap of the soli- 
ton continuum can be at incommensurate k, the lowest 
excited state almost always (actually except for very nar- 
row regions, whose existence may well be an artifact of 
our approximation) occurs at the commensurate values 
of k — or 7r/2 (see Fig. ||), which is consistent with 
the results of Refs. f^|ll|. However, the contribution from 
solitons can be observed in the dynamic structure fac- 
tor S(k,u). Unfortunately, the only existing dynamical 
DMRG studyE.il provides data only for one point in the 
dimerized phase (namely, for the point J s = j, J T = |, 
e = 0) where no incommensurability is expected. 

One can see that if e tends towards the multicritical 
points e = — 2 or e = +oo, then there are two transi- 
tions along the generalized symmetry line in the vicinity 
of the point J = | belonging to the dimer line: one at 
J < |, and another at J > |. This indicates that if the 
anisotropy exceeds a certain critical value, an additional 
transition boundary should appear in the (J T , J s ) plane. 
Note that in this case the variational phase boundaries 
lie close to the dimer line, so that this statement can be 
considered reliable. 



D. Comparison with numerical results 

We have studied the spectrum of low-lying excitations 
of the model (|J) by means of the exact diagonalization 
(Lanczos) method, using a chain consisting of 12 sites 
(i.e. having 12 spin and 12 pseudospin degrees of free- 
dom). 



Fig. |l| illustrates the behavior of the elementary exci- 
tations spectrum for a few choices of J s , J T , e, in com- 
parison with the variational results. For each value of 
the momentum k, the energies of four lowest-lying eigen- 
states in each of the subspaces with (S z ,t z ) = (0,0), 
(1,1), (1,0), and (0,1) were calculated. (The subspaces 
are denoted in Fig. |^ respectively with circles, diamonds, 
down and up triangles). One can see that our simple 
variational ansatz is in a qualitative agreement with the 
Lanczos data for the lowest excited states in the vicin- 
ity of the dimer line. It explains the main feature of 
the Lanczos data, namely, the rather flat dispersions ob- 
served nearly everywhere in the dimerized phase. One 
can also observe the tendency of solitons to have a spec- 
trum with the minimum at an incommensurate value of 
k when J s , J T decrease below certain value (see Fig. [|b) . 

Identification of different types of excitations can be 
done as follows: Free soliton-antisoliton pairs can have 
any of the four above-mentioned (S z ,t z ) combinations, 
whose energy would be exactly degenerate in the infinite 
system; that allows one to identify the soliton states in 
Fig. [|a-d as highly degenerate (in reality, nearly degener- 
ate) points where several different symbols are sitting on 
top of each other. S'-magnons appear as degenerate com- 
bination of circle and down triangle, r°-magnons should 
be circles, and r + -magnons must be up triangles. 

The agreement concerning the higher-lying excitations 
is considerably worse. Only for e = (see Fig. ^a,b) one 
can see the isolated magnon branch; for high anisotropics 
in the vicinity of the multicritical point (Fig. ^e,f) a large 
number of states appears immediately above the lowest- 
lying states appears, which may mean that actually many 
bound states split from the soliton continuum, and we are 
able to capture only the lowest-lying, most tightly bound 
states due to the short-range nature of the magnon ansatz 



IV. DOPED SPIN-ORBITAL MODEL: CHARGED 
SOLITONS AND THEIR BOUND STATES 

It is interesting to study the behavior of the spin- 
orbital model doped with holes (i.e., away from the quar- 
ter filling). In this case one has to consider the cor- 
responding t-J model, adding the appropriate hopping 
term to the Hamiltonian H(J S , J T ,e) determined by Eq. 

H = H{J S1 J T , e) + * X>U c *+i >CT r + • (16) 

lar 

Here q i(TT are the electron Fermi-operators at the site 
a and r are the corresponding spin and orbital quantum 
numbers, and the condition of no double occupancy is 
implicitly assumed. We also assume, in a similar fashion 
as in the course of derivation of (|I|), that hopping exists 
only between the nearest neighbors with the same type 
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of orbitals; the hopping amplitude t is chosen to be neg- 
ative for the sake of definiteness. One may also wish to 
add the nearest-neighbor Coulomb repulsion term 



H ^ H 



V J2 



nini+i 



J 



where n; = J2 aT 

Consider a single hole in this system. If one is on 
the dimer line J T — 3/4, J s = (3 + e)/4, where the 
checkerboard-type singlet state is exact, it is easy to show 
that the hole state \h;n) shown schematically in Fig. || 
gets hopped to the neighboring site without disturbing 
the structure of the ground state, so that the correspond- 
ing translational invariant state with momentum k is an 
exact eigenstate with the energy 



E lh {k) = 2e + 2<cos(/c), 



_ 3(3 +f) n7 x 
= -jg- • (17) 



Such a state connects two regions with different dimcr- 
ization pattern, so that it can be viewed as a "charged 
soliton" , or a hole bound on a soliton. One can straight- 
forwardly see that a hole state which does not disturb the 
dimer order is not able to move and has to decay into a 
"charged soliton" and an antisoliton. 

The two-hole problem can be also easily studied. The 
general two-hole state with a total momentum k can be 
written as 

*(fc)= e ik ^ +n ^ 2 f(n 2 -n 1 )\n 1 ,n 2 ), (18) 

and from the Schrodingcr equation one obtains the fol- 
lowing system of equations for the function f(r): 



2t cos - [/(r + 1) + f(r - 1)] = (E 2h - 4e )/(r), r > 2 , 



2tcos-/(2) = (E 2 h — 3e — V) /(l) . 



(19) 



Setting f(r) = Ae lqr + Be~ tqr , with the parameter q 
(having the meaning of the relative momentum) being 
real, one obtains the continuum of scattering states with 
the energy 

E 2h (k,q) =4e + 4£cos(fc/2)cos(< ? /2). (20) 

There is also a bound state solution with f(r) — foe~ Kr \ 
the parameter k > is determined by the equation 



(21) 



2tcos(fc/2) ' 

This solution describes a localized state with the energy 



E b 2hb (k) = 3e Q + V 



2i 2 



e -V 



(1+cosfc). (22) 



The spectrum of the 2-hole problem is schematically 
shown in Fig. @. One can see that for V < eo — 2\t\ 



the energy of this state lies below the lower boundary of 
the two- hole continuum in the whole range of k. When 
e < V < eo — 2\t\, the bound state exists only in a cer- 
tain interval [fco,7r] where ko is determined from (|^) by 
demanding that n is positive; ko tends to n as V tends 
to eo. For V > eo the energy of the localized state (|2^ ) 
lies above the upper boundary of the continuum, so that 
the states are "antibound." Thus, the "charged solitons," 
which are themselves bound states of a hole and a soliton, 
tend, in their turn, to form bound states, if the condition 
e > V is satisfied. 

One can observe that the general iV-hole problem is 
equivalent to that of the XXZ ferromagnet with J z = 
eo — V and J xy = 2t, in external magnetic field h = eo + V^ 
along the z axis, and thus can be treated by means of the 
Bethe ansatz. In fact, one can use the known result for 
the energy of TV-spin bound complexes to translate it into 
the corresponding formula for our problem: 



2 ,,21 1/2 



E N (k) = (e + V)N + {(e - V) 2 - U 2 } 



(23) 



x | tanh(iV//Vo) 



sinh(2/Vy/Vo) 



sin 2 (/fc/2) 



where A*o is defined by the relation 

cosh(2//V ) = (V- e )/2t. 



One can straightforwardly check that (|22j) is the partic- 
ular case of ( p3| ) for N = 2. For large N such a TV- hole 
complex corresponds to aj-classical soliton state in the 
XXZ ferromagnetic chainEj The expression (|2^) is valid 
for \V - e | > 2\t\. It is knownB that for \V - e | < 2\t\ 
only bound (respectively antibound) states with N = 2 
exist, and the energies of states with larger N fall into the 
continuum. For V — eo < —2\t\ the iV-hole bound com- 
plex is stable against the decay into smaller composite 
particles ,EMij i.e., 

E N (k) < min{E N _ M {k - k') + E M (k')} , 

and the opposite inequality (meaning respectively insta- 
bility) is fulfilled if V-e > 2\t\. Thus, for V-e < -2\t\ 
the ground state of the model (|T^) at finite doping rep- 
resents a condensate in which all holes areJaound in a 
single "particle" with zero total momentum;c3 in the in- 
termediate regime —2\t\ < V~eo < 2\t\ only the two-hole 
bound states are stable, and for V — cq > 2\t\ all bound 
states become unstable. 



V. SUMMARY 

We have studied, by means of the combination of vari- 
ational and exact results, as well as numerically, the in- 
terplay between different types of elementary excitations 
in the model of an anisotropic spin-orbital chain. We an- 
alyze the vicinity of the special "dimer line" in the phase 



G 



space of the model where the spontaneously dimerized 
ground state is known exactly. Solitons in the dimer 
order are shown to be the lowest excitations only in a 
narrow region of the phase diagram, and the phase tran- 
sitions are always governed by magnon-type excitations 
which can be viewed as soliton-antisoliton bound states. 
It is suggested that an additional phase boundary ap- 
pears when the anisotropy exceeds certain critical value. 

We have considered also a t-J version of the doped 
spin-orbital model. The elementary charge excitation is 
shown to be a "charged soliton" , or a soliton bound on 
a hole. We have studied the conditions under which the 
"charged solitons" tend to form bound states; the exact 
expression for the energy of the iV-hole bound complex, 
based on the known results from the Bethe Ansatz, is 
presented. 
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2n+l 



FIG. 1. The single soliton states |L)£ T , \R) a n T used in Eq. (|). The lower and the upper chains correspond to real spins 
(S) and orbital pseudospins (r), respectively. Thick solid lines denote singlet bonds. The soliton connects two degenerate 
spontaneously dimerized states which are exact ground states of the model (□) at J T — 3/4, J s — (3 + e)/4. 



"V > = 



• < 































• • • * = lt+t_>+ |t + t_>- (l+co)|t t > 

FIG. 2. Magnon states used in Eqs. dill), (|l^). Here w = for r-magnons (i.e., when the upper leg here is the r chain) and 
lo = e for S-magnons (the upper leg is the S chain). 



8 




FIG. 3. Variational results for the behavior of the gaps for different types of excitations, along the "generalized symmetry 
line" J s — J(l + e/3), J T = J, for a few values of the anisotropy e. Magnon gaps (shown by dashed lines) are straight lines. The 
gap of soliton-antisoliton pair excitations (shown by a solid line) is a straight line while the minimum of the soliton dispersion 
lies at k = 0, and it starts to deviate from a straight line for J below certain threshold value, when the minimum occurs at an 
incommensurate value of k. 
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FIG. 4. Behavior of the soliton and magnon excitations along the line J s = J T (1 + e/3) in the (J s , J T ) plane for fixed e. 
Solid lines indicate the boundaries of the soliton continuum, and isolated long-dashed lines correspond to the magnon-type 
excitations which can be viewed as soliton-antisoliton bound states. In (f), additional dashed lines near k = and k — it 
show the lower boundary of the magnon continuum formed by two r +1 magnons. Symbols denote the exact diagonalization 
(Lanczos) data for a system of the length 12 (rungs); circles, diamonds, up and down triangles correspond to the states with 
the quantum numbers (S z ,t z ) — (0,0), (1, 1), (0, 1), and (1,0), respectively. (For e = 0, up and down triangles are degenerate 
and thus for the sake of simplicity, only up triangles are displayed). 
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FIG. 5. A "charged soliton" state \h;n) corresponding to a hole localized on a soliton in the dimer order. 
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bound: V<e 
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FIG. 6. The spectrum of the two-hole problem, for fixed eo and t, at several different values of V, as discussed in the 
main text. Shaded region denotes the continuum of scattering states formed by a "charged soliton" and "charged antisoliton" ; 
isolated solid and dashed lines correspond to the bound and antibound states, respectively. 
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